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Abstract 

We compare the perturbatively calculated QCD potential to that obtained from lattice 
calculations in the theory without light quark flavours. We examine Et t(r) = 2m po \ e + 
^QCD( r ) by re-expressing it in the MS mass m = m MS (m MS ) and by choosing specific 
prescriptions for fixing the scale \i (dependent on r and m). By adjusting m so as to 
maximise the range of convergence, we show that perturbative and lattice calculations 
agree up to 3ro — 7.5 GeV -1 (tq is the Sommer scale) within the uncertainty of order 
A3 



QCD 
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1 Introduction 



For decades, the static QCD potential VqcdO"), formally defined from an expectation value 
of the Wilson loop, has been widely studied for the purpose of elucidating the nature of the 
interaction between heavy quark and antiquark. In modern language, a link to physical reality 
can be made naturally in the frame of potential-nonrelativistic QCD (potential-NRQCD) for- 
malism P3I21III]; in which Vqcd(^) is identified with the leading potential in expansion in 1/m 
of the heavy quarkonium system. Therefore, Vqcd(^) dictates, for instance, the bulk of the 
spectra of the bottomonium and charmonium states. 

Some time ago, there was a breakthrough that drastically improved the predictive power of 
perturbative QCD for the QCD potential and the heavy quarkonium spectrum: the perturbative 
predictions for these quantities became much more accurate. This was achieved by properly 
eliminating contributions from infrared (IR) degrees of freedom in the computations [HE] • The 
central quantity is the total energy of a static quark-antiquark pair, defined by the sum of 
the quark and antiquark pole masses and the QCD potential, E tot {r) = 2m polc + Vqcd(?")- We 
can achieve the decoupling of IR degrees of freedom (renormalon cancellation) at each order 
of the perturbative expansion by (1) re-expressing the quark pole mass in terms of a so-called 
short-distance mass, such as the MS mass, and (2) expanding m po i e and VqcdO") i n the same 
coupling constant.* As a result, the perturbative predictions become stable against a variation 
of the renormalisation scale n, and also the perturbative series show a much better convergence 
behaviour, as compared to those in the conventional computations. 

It was then natural to compare the perturbative QCD predictions with existing experimen- 
tal data or with other theoretical predictions which incorporate non-perturbative effects. The 
main aim of this program is to clarify the differences between the perturbative QCD predic- 
tions and the full QCD predictions, given the more accurate predictions of the former. The first 
comparison JB] was made for the bottomonium spectrum (and also for part of the charmonium 
spectrum) between the perturbative prediction and the experimental data. It was followed by a 
comparison [7j between the perturbative QCD prediction of the QCD potential and typical phe- 
nomenological potentials (used in phenomenological approaches to heavy quarkonium physics), 
and then by a comparison 8j of the QCD potential between perturbative QCD predictions 
and a lattice computation. More elaborated analyses on each of these comparisons followed 
subsequently 01101111111211121 ■ In all of these analyses, when IR contributions were appropri- 
ately eliminated, the perturbative QCD predictions turned out to agree with the experimental 
data/phenomenological potentials/lattice results within estimated perturbative uncertainties. 
Contrary to wide beliefs, there were no indications of large non-perturbative effects. Only much 
smaller non-perturbative contributions, which can be absorbed into perturbative uncertainties, 
appear to be compatible with these analyses. 

In this paper we are concerned with the third type of comparison: perturbative QCD predic- 
tion vs. lattice calculations. In the previous comparisons of this type, the leading renormalon 
uncertainty of the perturbative QCD potential was removed in various manners. In [Hj, the 
inter-quark force (improved by renormalisation group) was used instead of the QCD poten- 

*This is somewhat involved technically, since usually TO po i e and Vqcu(t) are expressed in terms of different 
coupling constants. 
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tiaF; in [TT] . the leading renormalon contribution (estimated by a sophisticated approximation) 
was subtracted from the QCD potential by hand; in [J2], the perturbative series was Borel- 
resummed, taking into account the leading renormalon pole appropriately We examine yet 
another method for removing the leading renormalon. Namely we examine the total energy 
Etot{r), as defined above, after re-expressing it in terms of the MS mass renormalised at the 
MS mass scale, m = m MS (m MS ). To achieve stable predictions over a wide range of r, we adopt 
the scale-fixing prescriptions of pHEI]. These prescriptions introduce the scales dependent on 
r and fn, /i — /x(r, fn), which are consistent with physical expectations. 

As stated, E tot (r) constitutes the leading part of the non-relativistic Hamiltonian of the 
heavy quarkonium system within the potential-NRQCD framework. We may expect that a 
direct comparison of E tot (r), which determines the bulk of the heavy quarkonium spectrum, 
would provide a clearer picture of the present status on the credibility of the theoretical predic- 
tions based on the potential-NRQCD framework, supplemented either by perturbative QCD 
computations or by lattice computations of the potentials; see e.g. |14 [ llf) | ll(i [ ITH] for analyses in 
this direction. Furthermore, the scale-fixing prescription for E tot (r) we adopt here is the only 
prescription which has been used in the perturbative QCD predictions for the level structure 
of the bottomonium states (including higher excited states) incorporating the renormalon can- 
cellation. We will see that indeed this prescription stabilizes the perturbative prediction up to 
large distances, and hence it is suited for predicting the energy levels of excited states of the 
heavy quarkonium systems. 

According to the renormalon argument, an uncertainty of the perturbative QCD prediction 
for E tot (r) grows rapidly at large distances as AQ CD r 2 jTT]. It is nevertheless important to 
predict E tot (r) perturbatively at large r for the following reasons. (1) The level spacings among 
the bottomonium spectrum have uncertainties smaller than the uncertainties of individual 
levels. This is because the errors of E tot (r) at different r are generally correlated. Indeed, 
the estimate of the error of E tot (r), by changing input parameters or scale-fixing prescriptions, 
is perfectly consistent with Aq CD r 2 [TUTU] ; on the other hand, the bottomonium level spacings 
vary less, because the individual levels vary in a correlated way. This is why the perturbative 
QCD predictions of the whole level structure of the bottomonium in jBUHlEl made sense. 
(2) Many physical quantities of heavy quarkonium states are sensitive to short-distance part of 
the potential. For instance, the fine splittings of the bottomonium excited states are sensitive 
to much shorter distance part of E tot (r) as compared to the individual levels. As a result, 
perturbative uncertainties of the fine splittings are much more suppressed (of order AQ CD /m 2 ) 
as compared to uncertainties of the individual levels which directly reflect uncertainties of 
Etot( r )- Predictability of E tot (r) up to large distances ensures that the wave functions of the 
excited states can be computed in the computation of the fine splittings [inj, although only 
the short-distance parts of the wave functions are relevant. The order Aq CD r 2 uncertainty of 
Etot{r) at large distances is just appropriate to ensure the theoretical uncertainties (Aq CD /m 2 ) 
of the fine splittings. These theoretical uncertainties, as well as the level of agreement with the 
experimental values, of the computed fine splittings (and the hyperfine splittings) turn out to 
be comparable to those of the recent lattice computations of these splittings; see |T3J[TB] for 

^See also [Jj for some theoretical discussion on the inter-quark force. 
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details. 

In our comparison of the QCD potential between perturbative QCD and lattice computa- 
tions, we benefit from considering a hypothetical world which contains no light quark flavours. 
It is then possible to use the lattice calculations of the QCD potential in the quenched ap- 
proximation. On the other hand, in the perturbative prediction for E tot (r), we have an ad- 
ditional parameter. Although naively the quark mass is simply a constant independent of r, 
due to our specific scale-fixing prescriptions, the value of m affects the r-dependence of E tot (r) 
non-trivially. For a heavy quarkonium system in this hypothetical world, there is no strong 
motivation to choose a specific value for m (as opposed to the studies [7J[10]). Therefore, in 
our analysis, we treat m as a controllable parameter for testing stability of the perturbative 
prediction. We will show that for those choices of m that give stable predictions, E tot {r) is 
independent of m up to deviations of the order of the expected theoretical uncertainty (after a 
suitable shift by an r-independent constant). By varying m to achieve optimum convergence 
for large r, we can obtain perturbative QCD predictions up to fairly long distances and compare 
them to the results of lattice QCD. 

The organisation of the paper is as follows: Sec. Insets our conventions and gives some details 
of our perturbative QCD calculation. Sec. El compares the lattice and perturbative QCD data. 
Conclusions are given in Sec. EJ We collect formulae related to the renormalisation-group 
evolution of the strong coupling constant in the Appendix. 



2 Conventions and framework 

We would like to compare the lattice data and the perturbative predictions corresponding to 
the same theoretical input. This will be carried out in the following manner. For each lattice 
data set we calculate the Sommer scale ro defined by 

2 dVQCD 



r 



dr 



1.65. (1) 



-ro 



Then the lattice data are expressed in units of r . The perturbative computations are expressed 
in terms of the strong coupling constant defined in the theory with n\ = active flavours. We 
convert all the results into units of ro using the relation between the Lambda parameter of the 
running coupling constant (in the MS scheme) and the Sommer scale pi)j: Aj^g = 0.602(48) Tq . 
We use the central value of this relation in the main part of our analysis; the effect of a variation 
of Ajyjg inside the error interval is discussed at the end of Sec. El All the predictions are compared 
in units of ro- Furthermore, in order to maintain physical intuition, we will also use physical 
units. Although there exists no rigid correspondence between the physical scales of the real 
world and of the hypothetical world, we follow the convention of the lattice calculations in 
the quenched approximation. The numerical value on the right -hand-side of Eq. has been 
chosen so that for phenomenological potentials ro ~ (400 MeV) -1 . Whenever we refer to values 
in units of MeV or GeV, we invoke this translation. 

The total energy of a static quark antiquark system is given by 

E tot (r) = 2m po i e + V QC d (r) . (2) 
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In perturbative QCD, the pole mass m pole is related to the MS mass m up to three loops by 
the relation 

_j 4a s (m) fa s (m)\ 2 ( a s (m)\ 3 . 
m pole = m | 1 + - — — + I ~ ^~ J <*i + M 2 r • (3 > 

The QCD potential up to 0(a s ) is given by 

4as(/i) 



VQCD(r) 



3 r 



+ 



An 
Air 



$l(At 2 + y) + 2(/3i + 2/3 a!)£ + a 2 



(4) 



where £ = log(yur) + je- 

Here and hereafter, we have set the number of light flavours, ni, to zero, i.e. we will be 
neglecting the effects of light quark loops. This corresponds to the quenched approximation 
in lattice QCD to which we want to compare the perturbative results. The running coupling 
Ois(ji) depends on ni through the coefficients of the beta function; the constants di,d 2 , a± and 
a 2 also get contributions from light quark loops and therefore depend on n\. For rii = 0, their 
values are (3 = 11, A = 102, d 1 m 13.443, d 2 » 190.39, a x = 31/3 and a 2 w 456.74. The 
analytical formulae can be found in |1U|.* 

After re-expressing as(m) in Eq. (jSJ) in terms of as{^) by using the running of as [see 
Eq. (JHJ) in Appendix], and dropping terms of 0(as{^) 4 ) and higher, we obtain the total energy 
Et ot (r;m,as(lJ>), A*) which does not suffer from the leading renormalon uncertainty. 

Due to the truncation of the perturbative series at finite order, E tot depends on the renor- 
malisation scale fi. Two scale fixing prescriptions have been introduced in [Jj: 

1. The scale \i = /ii(r) is fixed by demanding stability of E to t(r) against variation of the 
scale: 



d/j, 



0. 



(5) 



2. The scale /i = /i2(r) is fixed to the minimum of the absolute value of the last known term 



term] of E tot (r) 



d 



(6) 



Although these prescriptions are very different, it has been shown that where both prescriptions 
exist, the total energy is virtually identical for both prescriptions. As a general feature of £t t(r), 
the convergence of the perturbative series improves and the scale dependence decreases, if we 
choose larger /i for smaller distances and smaller \x for larger distances. Consequently, the range 

* These formulae have originally been computed in |2(JU21j . The mass relation (0 is re-expressed in terms of 
the coupling of the theory without heavy quarks. 
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0.022 
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2.8 
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1.140 


-0.102 


-0.029 
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1.006 


-0.005 





6.601 


3.0 


0.530 


1.179 


-0.133 


-0.021 


7.025 





576 


1.055 


-0.035 





7.019 


3.2 


0.543 


1.224 


-0.168 


-0.012 


7.444 





453 


1.680 


-0.682 





7.398 


3.4 


0.553 


1.277 


-0.208 


-0.004 


7.865 





507 


1.451 


-0.393 





7.859 


3.6 


0.559 


1.342 


-0.262 


0.006 


8.287 





577 


1.288 


-0.207 


0.005 


8.287 


3.8 


0.554 


1.439 


-0.352 


0.023 


8.711 





615 


1.268 


-0.173 


0.012 


8.708 



Table 1: Convergence properties of E to t for r = 2tq ~ 5 GeV . All numbers in GeV. 



of the perturbative calculation can be extended to much larger r with these prescriptions than 
what would be possible with a fixed, r-independent scale. 

The prescriptions for the renormalon cancellation and the scale fixing we adopt here follow 
(basically) those in [TjlTIHIT^]. We refer the reader to these papers for more detailed features of 
the perturbative predictions in these prescriptions. 

There are several methods to assess the reliability of the prediction at a given distance: e.g. 
one can compare the total energies as determined with the scales from both prescriptions, one 
can compare the sizes of the individual terms of the perturbative expansion of i^tot or one can 
study the scale dependence of E tot around the respective scale. We will use these methods in 
section 01 

The convergence properties of E tot (r) strongly depend on the mass parameter m. We 
illustrate this in Tabledfor r = 2r ~ 5 GeV -1 . E^ t denotes the 0{d l s ) term of the perturbative 
series of -Etot- It can be seen that the series converges nicely, especially for m ~ 3 GeV. In 
section El we will see that those values for m that provide best convergence also provide an 
optimal agreement with lattice results. We find that, to our surprise, the perturbative series 
converges for distances as large as 3ro — 7.5 GeV -1 . (Generally, the convergence behaviour of 
-E"tot(^) becomes worse for larger r.) We note that the values of the scales fii t2 stay much larger 
than 1/r. See [7j for discussion on this aspect. 

For large distances, -Etot depends on whether we use the analytical solution or numerical 
solution of the renormalisation-group equation for the running of the strong coupling constant 
as(fi) (see Appendix). We show the total energy for both types of running in Fig. [T]for two 
values of m. In accordance with our previous works [HE] we employ the numerical solution 
of the renormalisation-group equation below. Our results do not change qualitatively if we use 
the analytical running instead. 
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Figure 1: Comparison between analytical and numerical running. The solid lines correspond 
to m = 3 GeV ~ 7.5/r , the dashed ones to m = 4GeV ~ 10/r . 



To end this section, we briefly summarise the argument on the renormalons included in 
E tot (r). 

(1) After the O(Aqcd) renormalon is cancelled, the perturbative series of E tot (r) is estimated 
to have the following behaviour. The 0(a^ +1 ) term of the series expansion of E tot (r) 
behaves as const, x r 2 n\ (/?o«5/(6vr)) n n 35 ^ 2 for n ^> 1, where 5 = Pi/Pq- Because of the 
factorial n\, the series is only an asymptotic series, namely it diverges for large enough 
n. Hence, there is a limitation to the achievable accuracy of the perturbative prediction 
for E tot (r). It can be estimated from the size of the terms around the minimum, n m 
6n/ (p as)] this gives an uncertainty of order Aq CD r 2 . The behaviour of the series depends 
on the value of the expanding parameter as = as(fi), or equivalently, on the choice of 
the scale fi. The uncertainty, AQ CD r 2 , is independent of as(n) or fi, nonetheless. (For 
details, see (2211231 ■) 

(2) Based on the above argument, we may optimise convergence of the series by appropriately 
choosing the scale /i. In this case, even with the series expansion up to 0{a 3 s ) we may 
estimate the uncertainty of E tot (r) to be of order AQ CD r 2 from the size of truncated next- 
order term or from the scale dependence of E tot {r) around the optimised scale. Indeed, 
explicit numerical examinations of E tot {r) up to C(a|) support this argument 



(3) One may factorise the infrared part of E tot (r) using operator product expansion [2]. In this 
way, one may absorb the order AQ CD r 2 renormalon into a matrix element of an operator, 
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while defining a Wilson coefficient perturbatively that is free from the renormalon and 
dependent on the factorisation scale. More generally, one may separate E tot {r) into 
perturbative coefficients free from renormalons and non-perturbative parameters (matrix 
elements) including renormalons. (This factorisation is beyond the scope of the present 
paper.) 

3 Comparison of perturbative and lattice calculations 

For comparison with the perturbative predictions of the QCD potential as explained in the 
previous section, four different sets of lattice data calculated in the quenched approximation 
are used: those from [23 {J3 = 6.8), from [23] {j3 = 6.0), from |HI (0 = 6.0) and from [27] 
(6.57 < (3 < 6.92). All the lattice data have been corrected using the lattice Coulomb potential 
to match the continuum definition of the QCD potential at short distances. 

For comparison of the perturbative and the lattice data, we have to account for an r— 
independent additive constant that is not determined by lattice calculations. Since perturba- 
tive calculations are more reliable at small distances and lattice calculations are more reliable 
at larger distances, we adopt the following procedure: The different sets of lattice data are 
converted to physical units with the lattice spacing as given by the authors of the respective 
papers, or (where the lattice spacing was not explicitly derived) by fixing the Sommer scale 
with the phenomenological potential fit as performed by the authors of the respective papers. 
Then we adjust the sets of lattice data among each other to make them coincide at r = r 
by adding constants. Finally we shift both the perturbative and the lattice data so that they 
vanish at r = vq/4, where in the case of the lattice results the data from j2TJ is used. 

We see that the sets of lattice data PHl23[l2^ l2*7j corresponding to different values of (3 are 
located almost on the same curve (Fig.[2J). This shows that the dependence of the lattice results 
on the lattice spacing is negligibly small, i.e. discretization errors in the lattice calculations are 
negligible in our comparison. 

As described before, the perturbative calculation has two input parameters, these can be 
e.g. Ajyjg and m. The potential depends on the mass m (after shifting to -Ea ot (ro/4) = 0) only 
through the log(m/ /x)-terms in the relation between the pole mass and the MS mass. We find 
that (after shifting the curves to make them coincide at r /4), for small distances r < r /4 
the curves are identical, for larger distances the curves corresponding to different m start to 
differ from each other. The quality of convergence and the stability against scale changes varies 
strongly with m. 

In Fig.[2]we plot the QCD potential for various values of m between 1 and 6 GeV in steps of 
0.2 GeV. To ensure that only reasonably stable and reliable predictions are shown, the curves 
are drawn only in those points, where the energies as determined by the two different scale 
fixing prescriptions differ by less than 0.5/r . We find that the resulting curves span a band 
around the lattice data that increases in width with increasing r. The width of this band is 
consistent with the expected theoretical uncertainty due to the uncancelled next-to-leading 
renormalon [17J, ±|A 3 r 2 , with A = 300 MeV, indicated by the error bars in the figure. We find 
a very good agreement between the lattice results and the curves that show the largest range 
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Figure 2: Comparison between perturbative and lattice calculations of the QCD potential. 
The lines correspond to masses between 1 and 6 GeV in steps of 0.2 GeV (solid lines for 
integer masses, 1 GeV ~ 2.5/tq, the lines for m > 4 GeV are masked by the other lines). The 
points correspond to Bali/Schilling (o) Takahashi et al. (o) [2*5] . JLQCD (*) "2H] and 
Necco/Sommer (•) [22] . Error bars for the statistical errors of the lattice data are plotted 
where given by the authors, but they are generally smaller than the size of the symbols. Lines 
are plotted only when the total energies determined by the two prescriptions differ by less than 
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Figure 3: Comparison between perturbative and lattice calculations of the QCD potential. The 
lines correspond to masses between 1 and 6 GeV in steps of 0.2 GeV. Lines are plotted only 
when a +10% scale change makes the total energy vary by less than 20%. 



of convergence, but even for those choices of m where the prediction becomes unstable earlier, 
the agreement is still good. 

To show that the good agreement between the perturbative and the lattice calculations does 
not depend on a specific stability criterion, in Fig. El we show the same comparison as in Fig. El 
but this time we do not consider the difference between the energies as determined with the 
two different scale prescriptions, but the stability against scale change. We plot the curves only 
in those points where a scale change of +10% makes the total energy vary by less than ±20%. 

We would like to stress that we do not tune the mass parameter to achieve good agreement 
with the lattice results, but we vary it to find those values of m that give optimal convergence 
of the perturbative series. It can be seen in the figures that the curves for those values of m that 
have the largest range of convergence, the agreement with the lattice data is close to optimal. 

We now compare our results to those of [Hj. In that paper, a fixed, r-independent scale 
\i is used for the perturbative QCD potential. We find that our formalism almost exactly 
reproduces the curves of ^T] for large values of m (Fig. HJ. The explanation for this behaviour 
is the following: While in our formalism the scale is strongly dependent on r even for large 
masses (see Fig. |5J), the scale tends to rise with m. For m ~ 3 GeV the r-dependent scale 
varies around 1 GeV, for m ~ 10 GeV it varies around 3 GeV. Independent of fn, however, the 
scale dependence of E tot is strong for scales around 1 GeV and very weak for scales around 3 
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Figure 4: Comparison between our formalism and the 2-loop QCD potential according to the 
formalism of [TT]. The lowermost line shows the (m independent) potential for a constant 
fi = (0.15399 tq)" 1 , corresponding to the formalism of [IT]. The other lines show our results for 
masses from 1 to 12 GeV in steps of 1 GeV (solid lines for 1, 4, 7 and 10 GeV). 
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Figure 5: Scales as determined by the first prescription (0) for values of m between 1 and 12 
GeV. 




Figure 6: Dependence of E tot (r = 0.3 r ) on the scale fi for values of m between 1 and 15 GeV. 
Independent of m, the scale dependence is strong for fi < 1 GeV and weak for fi > 1 GeV. 
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r/r 

Figure 7: Dependence of the perturbative results on Aj^gr , the curves shown correspond to 
the centre and upper and lower bounds of the error interval given in ^H]- The middle curve 
corresponds to the one in Fig. El that has the error bars attached to it (m = 2 GeV). 

GeV (Fig. therefore choosing a large m in our formalism gives a result very close to the 
treatment of [TI]. We can also see that for these large masses our stability criteria indicate a 
range of convergence up to about r /2. Our analysis is therefore consistent with that of [11J, 
and the results of the latter are reproduced with our formalism by choosing large values for m. 

The perturbative predictions which are most stable at long distances in Figs.|2]and|B](m ~ 2- 
3 GeV) turn out to be less steep than the lattice data at r/ro ^ 0.5. Qualitatively, we expect 
that the larger the strong coupling constant the steeper the potential, because the interquark 
force becomes stronger This behaviour can be seen in Fig. where we have varied A^ro in 
the interval given in [HI, K m = (0.602 ±0.048) r \ The lower bound, centre and upper bound 
of this interval correspond to a!y ~°(M Z ) = 0.0801, 0.0811 and 0.08205, respectively. The larger 
value for Aj^g (if r ps 2.5 GeV -1 is fixed) results in a slightly steeper curve that reproduces 
the slope of the lattice data better than the central value. 

We would like to make two comments in this context: (i) We compared the perturbative 
QCD potential (including effects of light quark loops) with phenomenological potentials in |10j . 
There, the perturbative prediction with the input as(Mz) = 0.1181 (the present central value) 
turned out also to be slightly less steep than the phenomenological potentials at long dis- 
tances. As a result, a somewhat larger coupling as(Mz) = 0.1191-0.1201 was favoured for a 
better agreement with the phenomenological potentials, (ii) In [TT], the 0(ag) correction to 
the perturbative QCD potential (including the ultrasoft effects) was estimated and included 
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in a comparison with the lattice data. The estimated correction makes the potential some- 
what steeper. This is consistent with a naive expectation that such an effect is caused by an 
acceleration of the running of the coupling constant due to the 4-loop coefficient of the beta 
function. Thus, agreement of the perturbative potential with the lattice data may become even 
better than our present analysis when the full next order correction is calculated and included 
in the future. Furthermore, we have confirmed that the agreement between the perturbative 
and lattice data becomes even better if we use the 4-loop running of the coupling constant. 

4 Conclusions 

We have compared perturbative QCD and lattice QCD predictions for the QCD potential. We 
examined the perturbative QCD prediction for E tot (r) = 2m po i e + Vqcd(^) 5 taking specific pre- 
scriptions for fixing the renormalisation scale \x. We find that, by adjusting the mass parameter 
m, the perturbative prediction can be made stable up to distances r ~ 3tq. Whenever we 
obtain stable perturbative predictions, they agree with the lattice data within the uncertainty 
estimated from the residual renormalon of order A 3 r 2 . We emphasise that we do not tune m to 
fit the lattice data, but we tune m to achieve stability of the perturbative prediction, and then 
the agreement follows. Comparisons of perturbative QCD predictions and lattice data have 
been performed before e.g. in [HJII2], but only up to distances of 0.5 r and 0.9 r , respectively. 
If we take an optimal value of m, our prescriptions for the perturbative prediction of the QCD 
potential seem to give stable predictions to furthest distances among those examined so far. 

Our analysis provides a firmer ground to the analyses of jHlElEI]; which predicted the 
bottomonium energy levels up to the n = 3 states using (essentially) the same scale-fixing 
prescriptions as in the present analysis. § We note that the same conclusion could not be 
drawn directly from the previous comparisons |B1EJE] between the lattice and perturbative 
computations of the QCD potential, because the prescriptions adopted in those analyses have 
never been used in perturbative computations of the heavy quarkonium level structure including 
higher excited states. In the light of our present result, the scale-fixing prescription adopted 
here is optimal for stable predictions for the energy levels of excited states. Our result supports 
the estimates of theoretical uncertainties by the next-to-leading renormalon made in |51l5|ll3j. 

Whenever stable perturbative predictions are obtained, all the perturbative predictions 
with different prescriptions for subtracting the leading renormalon agree with one another and 
also with the lattice data, within the estimated uncertainty. In particular, our perturbative 
predictions for large m reproduce the G(ot%) perturbative prediction of [TT]. The fact that 
the different prescriptions have led to mutually consistent perturbative predictions of the QCD 
potential, endorses the consistency of the perturbative analyses. The comparisons between the 
perturbative and lattice data, together with other types of comparisons, provide evidence to 
the hypothesis that the perturbative prediction agrees with full QCD within the order A 3 r 2 
uncertainty. In particular, we consider the independent examinations corresponding to the 

§The value of m, which stabilizes the perturbative predictions for E to t(r) up to the furthest distance, lies 
between the bottom and charm quark masses. In this sense, we are in a lucky situation in the predictions of 
the charmonium and bottomonium spectra. 
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physical reality (jii = 4 with the non-zero charm quark mass ^U]) and to the hypothetical case 
(ni = 0) to be non-trivial cross checks with respect to the validity of the hypothesis. We may 
further make non-trivial tests of the hypothesis by varying the number of quark flavours and 
even the gauge group in comparing the perturbative and lattice calculations. 

It is quite surprising that the perturbative calculations turn out to give stable predictions 
up to such long distances [3r ~ (130 MeV) -1 ]. At the present stage the reason is unclear. 
Although we do not know a clear-cut criterion at which distance a perturbative QCD prediction 
should break down, possible arguments may be as follows. One point is that the relevant scale 
H for E tot (r) is not equal to 1/r but considerably larger. Another point that may be worth 
noting is that the system under consideration would be optimally suited for perturbative QCD 
computations. It is a colour singlet system having a localised spatial extent, so that we may 
expect the decoupling of IR degrees of freedom to be realised in a most natural way.' 



Appendix 

In this paper we expand the perturbative series of the QCD potential and the quark pole mass 
in the strong coupling constant defined in the MS scheme. The coupling constant obeys the 
renormalisation-group equation: 

/i 2 ^«s(/i)=/?(«5(/i)) = -«5(/i)E/5„(^) + • (7) 

We include the coefficients of the beta function up to 3 loops in our analysis", i.e. (3$ = 11, 
Px = 102, p 2 = 2857/2 (m = 0), and p n = for n > 3. 

In rewriting as(m) in terms of ag(/i) in the fixed-order expression of E tot (r), we use the 
perturbative solution of Eq. (|7J): 



a s (m) = at s (n) 



v m' V 7r / I 4 v m 7 8 v m 



7T 2 



(8) 



This relation is inserted to Eq. (jHJ) and the series expansion is truncated at order a s (/i) 3 . Then 
E tot (r) is given as a function of r, m, \i and &s{p)- 

The value of Qg(//) is determined by the renormalisation-group evolution from the input 
Lambda parameter defined in the MS scheme, A^; see e.g. ref. |2E1 for the definition of Aj^g. 
The renormalisation-group evolution of as{n) is calculated in two different ways in this pa- 
per. When we refer to "numerical running" , we solve the renormalisation-group equation (|7j) 
numerically. In this case, A^ = 0.65573... x /^Landau ( n i = 0), where /iLandau is the position of 
the infrared singularity (Landau singularity) of the running coupling constant as{f^)- On the 

tit may be contrasted with e.g. perturbative QCD calculations of parton scattering amplitudes, where spa- 
tially separated coloured partons exist as asymptotic states. 

"Although the 4-loop coefficient is available, we consider the use of the 3-loop beta function to be more 
consistent with the analysis of the fixed-order perturbative series up to O(ajg). 
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other hand, when we refer to "analytical running" , we use an approximate analytic solution of 
the renormalisation-group equation: 



vr ~ (3 L ((3 L) 2 ' (J3qL) 



+ 



2 



4 



(log 2 L-logL-l)+/3 2 



(9) 



where L = log(yU 2 /A|^) and terms of 0(1/ L 4 ) have been neglected. 
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